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ABSTRACT 

We show that stars with transiting planets for which the stellar obliquity is large are preferentially hot 
(Tgff > 6250 K). This could explain why small obliquities were observed in the earliest measurements, which 
focused on relatively cool stars drawn from Doppler surveys, as opposed to hotter stars that emerged later from 
transit surveys. The observed trend could be due to differences in planet formation and migration around stars 
of varying mass. Alternatively, we speculate that hot- Jupiter systems begin with a wide range of obliquities, 
but the photospheres of cool stars realign with the orbits due to tidal dissipation in their convective zones, while 
hot stars cannot realign because of their thinner convective zones. This in turn would suggest that hot Jupiters 
originate from few-body gravitational dynamics, and that disk migration plays at most a supporting role. 

Subject headings: planetary systems — planets and satellites: formation — planet-star interactions — stars: 
rotation 



1. INTRODUCTION 

There are now 28 cases of stars with transiting planets 
for which the stellar obliquity — or more precisely its sky 
projection — has been measured via the Rossiter-McLaughlin 
effect. The history of these measurements is perplex- 
ing. Starting with the pioneering measurement of Queloz et 
al. (2000), for 8 years a case was gradually building that the 
orbits of hot Jupiters are always well-aligned with the rotation 
of their parent stars. Then in a sudden reversal, several mis- 
aligned systems were found, with the first sighting by Hebrard 
et al. (2008) and the most recent spate of discoveries by Tri- 
aud etal. (2010). 

In this Letter we point out that the misaligned systems are 
preferentially those with the hottest photospheres. In § 2 we 
discuss the sample, and in § 3 we display the patterns involv- 
ing the order in which the measurements were made, the stel- 
lar effective temperature, and the stellar obliquity. In § 4 we 
speculate on the meaning of the patterns, and in § 5 we sum- 
marize the results and their implications for theories of the 
origin of hot Jupiters. 

2. THE SAMPLE 

We focused on those systems for which the projected spin- 
orbit angle. A, was measured with a Icr precision of 10° or 
better The less precise cases are not as helpful because we 
cannot tell definitively whether the system is aligned or mis- 
aligned, and because the large uncertainties are usually asso- 
ciated with strong systematic effects. 

We omitted Kepler-^ (Jenkins et al. 2010) from considera- 
tion even though the quoted uncertainty is smaller than 10°, 
because no data were gathered immediately before or after the 
transit, precluding tests for a systematic velocity offset on the 
transit night. Such offsets are possible, or even probable, for 
stars as faint as Kepler-% observed in bright moonlight (see, 
e.g., Tripathi et al. 2010). When we reanalyzed the Kepler-% 
data allowing for such an offset, the result was A = 20° ±20°. 
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Table 1 summarizes the resulting sample of 19 systems, 
along with the properties of the 9 omitted systems, for com- 
pleteness. For simplicity we refer to the planets as "hot 
Jupiters" because they are all giant planets with short peri- 
astron distances, although it should be remembered that they 
span a wide range of masses (0.36-11. 8 Mjup) and orbital pe- 
riods (1.3-11 1 d). 

3. THE PATTERN 

The top panel of Fig. 1 shows A as a function of the date of 
the earliest reported observation of the Rossiter-McLaughlin 
effect. The trend of low values (good alignment) for the first 
few years is evident, as is the "spike" of high values (mis- 
alignment) in the most recent years. This plot also suggests 
that the systems initially discovered in radial-velocity (RV) 
surveys are systematically more well-aligned than those sys- 
tems discovered in transit surveys. 

The reason for this pattern may be that the earlier measure- 
ments focused on cooler and less massive stars. The bottom 
panel of Fig. 1 shows that the average effective temperature 
(Teff) of the host stars has risen with time. Only in 2008 did 
investigators begin examining stars with T^s > 6250 K, and all 
of those systems were identified in transit surveys as opposed 
to RV surveys. 

We cannot give a deterministic explanation for this trend, as 
it depends not only on the selection functions for the various 
surveys but also sociological factors affecting the allocation 
of telescope time. However it seems probable that cooler stars 
were examined earlier because they allow for better RV pre- 
cision, and therefore greater ease of confirming the existence 
of planets. Indeed, most RV surveys exclude early-type stars 
altogether. In contrast, transit surveys have nearly magnitude- 
limited samples that include hot and luminous stars. These 
factors may explain why planets around hot stars were only 
found in transit surveys, and why they emerged relatively late 
from those surveys. 

Fig. 2 shows A as a function of Teff. Most of the misaligned 
systems are around the hottest stars in the sample. The tran- 
sition from aligned to misaligned occurs around Tgff = 6250 K 
(spectral type F8), which for the rest of this Letter we take 
to be the boundary between "cool" and "hot" stars. We will 
also use the term "misaligned" to mean |A| > 10° with >3(j 
confidence. 
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Fig. 1 . — A brief history of Rossiter-McLaughlin measurements. Top. — The projected obliquity (spin-orbit angle) is plotted as a function of the earliest 
date of observation reported in the literature. Blue squai'es indicate systems discovered by RV surveys, while red circles indicate systems found in photometric 
transit surveys. Bottom. — The stellar effective temperature of the same systems. Systems with T^ff higher than 6250 K have been discovered by transit surveys, 
and began to be examined for the RM effect in 2008. 



Among the cool stars, two out of 1 1 (18%) are misaligned, 
while among the hot stars, 6 out of 8 (75%) are misaligned. 
Another way to describe the pattern is to enumerate excep- 
tions to the rule that only hot stars are misaligned. There are 
two types of exceptions: "strong" exceptions in which a cool 
star is misaligned, and "weak" exceptions in which a hot star 
is apparently well-aligned. Weak exceptions are not as se- 
rious because only the sky-projected obliquity is measured, 
and consequently a low value of A could be observed for a 
misaligned system. Out of 19 systems, there are two strong 
exceptions and two weak exceptions. 

Many seemingly compelling trends of this kind turn out to 
be spurious. The best way to make progress is to gather more 
data. The prediction that misaligned systems are preferen- 
tially around hot stars will be tested in the near future, and a 
primary purpose of this Letter is to enunciate the prediction in 
advance of forthcoming observations. 

It is also important to consider selection effects. From a 
transit surveyor's perspective, the most important difference 
between a well-aligned star and a misaligned star of the same 
spectral type is that the well-aligned star has a larger vsin/ 
(sky -projected rotation rate). A large vsin/ implies broader 
spectral lines and poorer Doppler precision, inhibiting planet 
discovery. Therefore there is a potential bias against discov- 
ering well-aligned systems. We must ask whether there could 
exist a large population of well-aligned systems around hot 
stars that has been missed by cuiTent surveys. 

For this question, the RV surveys are iiTelevant because they 
exclude all hot stars regardless of vsin/. As for the transit 
surveys, to assess the bias we must know how transit candi- 



dates are identified and followed up. Latham et al. (2009) pro- 
vided a complete inventory of transit candidates and follow- 
up observations, which we take to be representative. They 
chose 28 transit candidates for spectroscopic follow-up, with- 
out regard to spectral type. Of those, 4 were not pursued fur- 
ther once it was found that vsin / > 50 km s"'. The other 
24 cases were observed assiduously until a hot Jupiter was 
confirmed (2 cases) or ruled out (22 cases). Hence, any bias 
against well-aligned systems is probably only for stars with 
vsin / ^ 50 km s~' . Such rapid rotators typically have spectral 
types <F3 and Tgff > 6700 K, whereas our sample ranges from 
5040 to 6700 K, with one exception^ Hence it seems unlikely 
that this bias is completely responsible for the observed X-T^{f 
relation, although it may play some role. 

Another caveat is that we cannot tell whether the relevant 
parameter is really Tgff or some other correlated variable, such 
as stellar mass. The transition temperature of 6250 K coiTe- 
sponds to approximately 1.3 Mq for solar- me tallicity main- 
sequence stars. 

4. POSSIBLE EXPLANATIONS 

Despite these caveats it is impossible to resist speculating 
on the reasons why hot stars with hot Jupiters have high obliq- 
uities. We restrict ourselves to an airing of issues and a toy 

^ The exception, WASP-33 (7200 K), proves the rule. The planet was dis- 
covered despite the star's rapid rotation (I'sin; = 86 km s"') by exploiting 
the RM effect (Collier Cameron et al. 2010) and not by the usual proce- 
dure of measuring the spectroscopic orbit. Therefore, even if the other hot 
systems were selected in a manner biased against well-aligned systems, we 
would expect WASP-33 to be more representative of the true obliquities of 
hot stars — and it is misaligned. 
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Fig. 2. — Misaligned systems liave liotter stars. Top. — The projected obliquity is plotted against the effective temperature of the host star. A transition 
from mainly-aligned to mainly-misaligned seems to occur at T^fi 6250 K. The two strongest exceptions are labeled. Symbol colors and shapes have the same 
meaning as in Fig. 1. Bottom. — The mass of the convective zone of a main-sequence star as a function of Tiff, from Pinsonneault, DePoy, & Coffee (2001). It is 
suggestive that 6250 K is approximately the temperature at which the mass of the convective zone has bottomed out. 



model illustrating a speculative hypothesis, leaving detailed 
investigations for future work. 

One possibility is that there are two pathways for produc- 
ing hot Jupiters, one of which is specific to low-mass stars 
and yields low obliquities, while the other occurs mainly for 
massive stars and produces a broad range of obliquities. The 
low-obliquity mechanism could be insptral due to tidal inter- 
actions with the protoplanetary disk (Lin, Bodenheimer, & 
Richardson 1996). The high-obliquity mechanism could be 
some combination of planet-planet scattering (Chatterjee et 
al. 2008) and Kozai cycles (Fabrycky & Tremaine 2007). It 
is not obvious why these mechanisms would have a strong 
dependence on stellar mass or temperature, although it is in- 
teresting that 1.3 Mq is approximately the same stellar mass 
above which giant planets are found to have larger masses, 
wider orbits, and a higher rate of occurrence (Bowler et 
al. 2010). Perhaps more massive stars are more likely to form 
systems of massive planets in unstable configurations, leading 
to an enhanced rate of gravitational scattering in comparison 
to cooler stars. 

Another possibility is suggested by the sharpness of the 
transition from aligned to misaligned, and its location at 
Teff ~ 6250 K. For main-sequence stars, this is approximately 
the temperature above which the mass in the outer convective 
zone (Mcz) becomes inconsequential. The decline in Mcz is 
illustrated in the bottom panel of Fig. 2, based on the relation 
presented by Pinsonneault, DePoy, & Coffee (2001). Between 
spectral types GO and F5 (5940 and 6650 K), increases by 
a factor of 1.3, and Mcz decreases by a factor of 120. 

Convective zones are important for the production of mag- 
netic fields and for tidal dissipation. Magnetic fields may be 
relevant by setting the inner radius of the protoplanetary gas 
disk, where accreting material is captured onto field lines, or 
by allowing the star to spin down through magnetic braking. 



The possible relevance of tidal dissipation is even more obvi- 
ous, as it would tend to realign the star with the orbit. 

Pursuing this latter point, we hypothesize that there is a sin- 
gle mechanism for producing hot Jupiters, and this mecha- 
nism yields a broad range of obliquities. For the cool stars, 
tidal dissipation damps the obliquity within a few Gyr, while 
for the hot stars, dissipation is ineffective. Therefore we ob- 
serve hot Jupiters to be well-aligned around cool stars, and 
misaligned around hot stars. 

It has been argued previously that tidal dissipation is 
too slow to affect the stellar spin state (see, e.g., Winn et 
al. 2005), but these arguments should now be reconsidered. 
The timescales for tidal dissipation are not understood from 
first principles and are poorly constrained by observations. 
Another objection is that obliquity damping should be accom- 
panied by spin-orbit synchronization, which is not observed. 
However, cool stars spin down due to magnetic braking. Thus, 
even if tides do synchronize the rotation and orbital periods 
while damping the obliquity, magnetic braking could subse- 
quently slow the rotation to the observed values. A third ob- 
jection, and the hardest to overcome, is that obliquity damp- 
ing is accompanied by orbital decay, threatening the planet 
with engulfment (Levrard, Winisdoerffer, & Chabrier 2009, 
Barker & Ogilvie 2009). The planet must surrender all its an- 
gular momentum in order to reorient the star, because of the 
star's large moment of inertia. 

We are thereby led to explore a scenario in which the star's 
moment of inertia is drastically reduced. We suppose that 
only the convective zone is dissipatively torqued by the planet, 
and that the radiative zone is weakly coupled to the convec- 
tive zone and to the planet. Without the burden of the massive 
radiative interior, the convective zone — and thus the observ- 
able photosphere — can align with the planetary orbit without 
drawing in the planet. Likewise, the magnetic braking torque 



4 



Winn, Fabry cky, Albrecht, Johnson 2010 



would be even more effective in slowing the surface rotation 
speed and preventing spin-orbit synchronization. 

Core-envelope decoupling has been discussed in the context 
of young stars (see, e.g., Irwin & Bouvier 2009), but here 
we would need decoupling to persist for a sizable fraction of 
the main-sequence lifetime of a cool star A problem with 
this notion is that the Sun's convective and radiative zones 
appear to be well-coupled (Howe 2009). However, this may 
not have always been so, and it was not a foregone conclusion 
theoretically (see, e.g., Pinsonneault et al. 1989). The most 
plausible solar coupling mechanisms, magnetic linkage and 
internal gravity waves, may be absent or may act on longer 
timescales for stars with hot Jupiters. 

To investigate the effects of core-envelope decoupling we 
used the equations of Eggleton & Kiseleva-Eggleton (2001) 
to follow a circular orbit of a hot Jupiter around a 1 Mq star, 
with initial periods Porb = 3 d and Piot = 10 d. Based on the 
stellar evolution code EZ-Wet0 we take the convective zone 
to have mass 0.015 Mq, moment of inertia 0.0066 MqR^, 
and apsidal motion constant 9 x 10""^. We chose a tidal dissi- 
pation factor 2* = 6 X 10^, which is consistent with the cur- 
rent population of hot Jupiters, although the large uncertainty 
in 2* causes a correspondingly large uncertainty in all of the 
timescales reported here. We do not model the dissipative 
shear or the non-dissipative oblateness coupling between the 
convective zone and the radiative interior The magnetic brak- 
ing torque was modeled with an extra term in the equations of 
motion: 



—J— =-amb"*"*- 
at 



(1) 



For the braking coefficient a^b we used 1 .66 x 10"'^ yr, based 
on a scaling of the Barker & Ogilvie (2009) results according 
to the moment of inertia. 

Fig. 3 shows the time history of fiot, i^orb, and the stellar 
obliquity ip, assuming an initial value of 60 deg. (Similar 
results were obtained from an initially retrograde condition.) 
Three lines are plotted, corresponding to planet masses of 3, 
1, or 1/3 Mjup. For Jupiter-mass planets, the stellar obliquity 
damps before the planet is consumed. Magnetic braking pre- 
vents synchonization of the convective zone with the orbit, in 
agreement with observations. However, this model also im- 
plies that orbits decay within main-sequence lifetimes, and 
that close-in massive planets should be rarer around cool stars 
than hot stars, due to their more rapid orbital decay. 

Another prediction is that the planets exerting the weak- 
est tidal torques should be seen as "strong exceptions": mis- 
aligned planets around cool stars. To compute the obliquity- 
damping component of the tidal torque we averaged together 
the last terms of Eqns. (10) and (1 1) of Eggleton & Kiseleva- 
Eggleton (2001), giving a decay timescale proportional to 



(\-e 



2)9/2 



l + 3e2 + (3/8)e4' 



(2) 



where Mp is the planet mass, a is the orbital distance, is 
the stellar radius, and e is the orbital eccentricity. By this 
standard, the 3 systems with the longest timescales for obliq- 
uity damping are HD 80606, HD 17156 and WASP-8. Thus, 
in our theory it is appropriate that HD 80606 and WASP-8 are 
strong exceptions. 

' http://www.astro.wisc.edu/~townsend/static.php?ref=ez-web 
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Fig. 3. — Toy model in which obliquity damping precedes orbital de- 
cay. See the text for details. Each panel shows curves con'espondiiig to 
planetary masses of 3, 1, and 1/3 Mjup. Top. — The obliquity ij} of the observ- 
able photosphere damps quickly, because the outer convective zone has little 
mass. Bottom. — The star's rotation period, Prot,* (solid Hues), is kept slow 
by the action of magnetic braking. The planet's orbital period, Poib (dashed 
lines), ultimately shrinks to zero (and the planet is destroyed), but on a longer 
timescale than obliquity damping. 



5. DISCUSSION 

The finding that hot stars with hot Jupiters tend to have 
high obliquities is not the only pattern that has been described 
in the Rossiter-McLaughlin data. Johnson et al. (2009) and 
Hebrard et al. (2010) found that the first 3 known misaligned 
systems all involved relatively massive planets on eccentric 
orbits. Since then, several exceptions have been discovered, 
such as WASP-15 and WASP-17 (Triaud et al. 2010). 

The A-Feff relation may be a sign that the mechanisms that 
produce hot Jupiters depend strongly on stellar mass. We have 
also explored a theory in which hot Jupiters are emplaced with 
a wide range of obliquities around all stars, but the cool stars 
tidally realign with the planetary orbits. The main difficulty 
with any theory of tidal realignment is avoiding orbital de- 
cay. Core-envelope decoupling could postpone orbital decay 
until after alignment is achieved, although this scenario is ad- 
mittedly speculative. One implication would be that close-in 
massive planets should be rarer around cool stars. Another 
implication would be that attempts to compare the ensemble 
results for A and the predictions of migration theories, such 
as those of Fabrycky & Winn (2009) and Triaud et al. (2010), 
should consider only hot stars, because cool stars may have 
been affected by subsequent tidal evolution. 

Finally, we interpret the results, as did Triaud et al. (2010), 
as a blow against the theory of disk migration, which would 
yield low obliquities as a general rule. Disk migration prob- 
ably does play a role in sculpting exoplanetary orbits, and 
convergent migration of multiple planets may occasionally 
produce tilted orbits (Yu & Tremaine 2001). But if obliq- 
uity truly depends on the present-day convective zone of the 
host star, then hot Jupiters likely arrived after the pre-main se- 
quence convective phase ceased, tens of Myr after disk disper- 
sal. Few -body gravitational dynamics (scattering or Kozai cy- 
cles) followed by tidal dissipation in the planet is compatible 
with this timescale, and it naturally produce misalignments, 
so this mechanism might account for most or all hot Jupiters. 
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TABLE 1 

A BRIEF HISTORY OF ROSSITER-MCLAUGHLIN MEASUREMENTS 



TT 

Name 


Earliest observation date 


— 

Type or survey 






Tcff [K] 




A [deg] 


References 


HD 209458 


2000 Jul 29 


RV 




6070 ± 50 




-4.4 ± 1.4 


1,2 


HD 149026* 


2005 Jun 26 


RV 




6160±50 




-12±15 


3 


HD 189733 


2005 Aug 21 


RV 




5040 ± 50 




-0.85+0-2? 


4,5 


TrES-1* 


2006 Jun 21 


Transit 




5230 ± 50 




30±2I 


6 


TrES-2* 


2007 Apr 26 


Transit 




5850 ±50 




-9 ±12 


7 


HAT-P-2* 


2007 Jun 06 


Transit 




6290 ±60 




2+12-2 

'-'■^-12.5 


8,9 


HAT-P-1 


2007 Jul 06 


Transit 




6000 ± 120 




3.7±2.1 


10 


Corot-2 


2007 Jul 16 


Transit 




5600 ± 120 




-7.2 ±4.5 


11 


lrJ:io-4 


if\r\n Till 11 
zUU/ Jul I J 


Transit 




Aonn -1- n< 




—D.J it 4. / 


1 

Iz 




2001 Nov 1 2 


RV 




6080 + 56 




10.0±5.1 


13 14 15 16 


XO-3 


2008 Jan 28 


Transit 




6430 ±50 




-37.3 ±3.7 


17,18 


Corot-l* 


2008 Feb 27 


Transit 




6000 ±150 




77 ±11 


19 


HAT-P-7 


2008 May 30 


Transit 




6350 ±80 




182.5 ±9.4 


20,21 


WASP-3 


2008 Jun 18 


Transit 




6400 ±100 




3 3+2-5 


22,23 


WASP-18 


2008 Aug 21 


Transit 




6400 ±100 




-5.0+1] 
-37 6+22-3 

■^'■"-IQ.O 


24 


Corot-3* 


2008 Aug 26 


Transit 




6700 ± 140 




25 


WAor-o 


ZUU5 UCt U4 


Transit 




^Ann -1- cn 
joUU it oU 




1 14 -i-i-a.g 

-114.2_4g 


ZO 


WASP-4 


2008 Oct 08 


Transit 




5500 ±150 




4-1-34 
—♦3 


24 




ZUUo UCL Uo 


Transit 








11+18 


Z / 


WASP-2* 


2008 Oct 15 


Transit 




5200 ± 200 






24 


WASP-5 


2008 Oct 16 


Transit 




5700 ± 150 




1 2 4+^-2 


24 


WASP- 15 


2009 Apr 27 


Transit 




6300 ± 100 




-139.6+« 


24 


WASP- 17 


2009 May 22 


Transit 




6600 ± 100 




-147.3+1= 
42±8' 


24,28 


HD 80606 


2009 Feb 13 


RV 




5570 ±44 




29,30,31,32 


WASP-14 


2009 Jun 17 


Transit 




6500 ± 100 




-33.1 ±7.4 


33 


Kepler-8* 


2009 Oct 29 


Transit 




6200 ± 150 




-26.9 ±4.6 


34 


WASP-33 


2009 Dec 08 


Transit 




7400 ±200 




-107.7 ±1.6 


35 


HAT-P-13 


2009 Dec 27 


Transit 




5640 ±90 




-0.9 ±8.5 


36 


Note. — References 


: (1) Winn el al. (20Il.Sj. i2j Qiieki/ cl al. (2000). (.)) WcKcl al, (2007). (4) Tri 


and cl al. (2009). (5) Winn 


el al. (2006). (6) 




iWiimelal 


(2008). (8) Locillcl el al. (2008). (9) Winn cl al. (2007) (10) Johnson c 


al. (200S),(ll)BouehyelLil. (200S), (12} Nai ila el al. 12010). ( L^) Narila el al. {2009a). (14) Bai bieri el 


:al.(2009).(l.S)Cochiancl 


al, (2008). (16) Narila el al, (2008). 1 17 


iWiiiiicl al 


(2009a). ( ]8) llcbrard cl al. (2008), (19) Poni cl al. (2010), (20) Winn c 


al. (2009e). (2l)Narilac 


;tal. (2009b(, (22) Tripalhi el al. (2010). (2.3) Simpson el al. (2010). (24) Triaiid cl al. (20 10). l2.S) Ti iaLid 


el al, (2009). (26) Qneloz el al. (2010). i 


27iGilloiie 


lal.)2009). (28) Anderson el 


al. (2010). (29) llcbrard el al. (2010), (3C 


Moulouelal. (2009). (31 


1) Winn cl al. (2009b). (32) Pom el al. (2010). (33) Johnson el al. (2009). (34) J 


enkinsclal.(2010).(3.S)C( 


illier Cameron el 


al. (2010). (36) Winn e 


•I al, (21)11)1. 


Where more than one referen 


cc is given, ihc L|uolcd value for A is lake 


from the fir^l reference ii 


a the list. Some authors use a different coordinate system and report /3 = -A ; 


for this table we have convt 


;rted all results to 


A. For WASP-33 the tabulated vali 


ae and error bar for A represent the mean and standard deviation of the 


independently derived values given by Collier Cameron et al. (2010). Staned systems (*) wei« omitted from the sample discussed i 


n H 2-3. 











